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Abstract: Many search strategies for the Standard Model Higgs boson apply specific 
selection criteria on hadronic jets observed in association with the Higgs boson decay 
products, either in the form of a jet veto, or by defining event samples according to jet 
multiplicity. To improve the theoretical description of Higgs-boson-plus-jet production (and 
the closely related Higgs boson transverse momentum distribution), we derive the two- loop 
QCD corrections to the helicity amplitudes for the processes H — > ggg and H — > qqg in 
an effective theory with infinite top quark mass. The helicity amplitudes are extracted 
from the coefficients appearing in the general tensorial structure for each process. The 
coefficients are derived from the Feynman graph amplitudes by means of projectors within 
the conventional dimensional regularization scheme. The infrared pole structure of our 
result agrees with the expectation from infrared factorization and the finite parts of the 
amplitudes are expressed in terms of one- and two-dimensional harmonic polylogarithms. 
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1. Introduction 

Within the Standard Model of particle physics, the Higgs boson is the only particle re- 
maining to be discovered. The Higgs boson is crucial for electroweak symmetry breaking, 
the mechanism that explains the generation of the masses of the fermions and the weak 
gauge bosons. While the vacuum expectation value of the Higgs field is directly related 
to the Fermi constant, its mass remains a free parameter that can be constrained but not 
predicted by the theory. 

The direct detection of the Higgs boson at LEP and the Tevatron has been a very 
challenging task over the past two decades [1,2]. The LEP experiments [1] excluded Higgs 
boson masses below Mh ~ 114 GeV, while the Tevatron excluded Higgs masses in a 
narrow window around the VF-pair threshold Mh ~ 2M\y. With the start of the Large 
Hadron Collider (LHC) at CERN, the focus of the Higgs boson searches has moved to the 
ATLAS and CMS experiments, where the search is based on different decay channels. For 
Mh > 135 GeV the decay into two weak gauge bosons is most prominent, while for lower 
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Higgs boson masses the search at the LHC is much more challenging, since the dominant 
decay modes are overwhelmed by large Standard Model backgrounds. For example, the 
search for light Higgs bosons with Mh ^5 130 GeV is based on the rare decay H — > 77, which 
has a branching ratio of O (lO -3 ) , thus requiring larger integrated luminosity. For light 
Higgs boson masses, the dominant production process at the LHC is gluon fusion. Based 
on the first 0(5 fb -1 ) of proton data taken in 2011, ATLAS [3] and CMS [4] are now able 
to narrow down the allowed mass range for the Standard Model Higgs boson considerably 
by essentially excluding the Higgs bosons in the range 0(130 GeV) > Mh > 0(600 GeV), 
while observing an excess of Higgs boson candidate events around Mh = 125 GeV. 

At leading-order (LO), the Higgs coupling to the two gluons is mediated through a 
quark loop. Since the Higgs coupling to the quarks is proportional to the quark masses, 
the dominant contribution is generated from the top quark [5]. The next-to-leading-order 
(NLO) corrections [6] to this process have been calculated and turn out to be very large 
(>60%). In the heavy top quark limit, M t — > 00, the Hgg coupling becomes independent 
of M t . One can therefore integrate out the top mass (M t ) and formulate an effective 
Langrangian C e ff for the Hgg coupling [7]. This technique is valid for Mh < 2M t and 
reduces the loops that need to be calculated by one. In this limit, the inclusive Higgs 
boson production cross section has been computed at NLO [8] and at next-to-next-to- 
leading-order (NNLO) [9], indicating a stabilization of the perturbative prediction at this 
order. 

Experimental searches of the Higgs boson apply final state cuts to improve the sig- 
nificance of a potential signal over Standard Model background processes. To implement 
these cuts in the theoretical description, fully exclusive calculations, which keep track of 
the kinematical information of all final state particles (Higgs decay products and QCD 
radiation) are mandatory. In the heavy top quark limit, the NNLO corrections to Higgs 
production via gluon fusion have been computed fully exclusively, including the Higgs de- 
cay to two photons or two weak gauge bosons by two independent groups [10, 11]. These 
calculations are in the form of flexible parton-level event generators, which can properly 
account for the final state restrictions used in the experimental studies. 

An important final state discriminator is the number of jets observed in addition to 
the potential Higgs boson decay products, and the Higgs signal can often be enhanced 
by applying jet vetos [12, 13]. In many searches, it is however expected that the H + Oj 
and H + lj samples contribute roughly equally to the sensitivity. In the above-mentioned 
NNLO calculations, the H + lj final states are included to NLO [14], and the H + 2j 
final states to LO. NLO corrections to H + 2j-production have been derived recently [15]. 
The correlation between samples of different jet multiplicity has recently been matter of 
quite some debate [16,17], and an improved theoretical description of H + 1 j-production to 
NNLO accuracy is essential in order to have the same theoretical accuracy for the H + lj 
contribution as for the H + Oj contribution. 

In the heavy top quark limit, a full NNLO QCD calculation of H + lj production 
requires the computation of the matrix elements of three contributions: 

(a) the tree level H — > 5 partons amplitudes, 
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(b) the one-loop corrections to the H — > 4 partons amplitudes, 

(c) the two-loop corrections to the H — > ggg and H — > qqg matrix elements. 

The tree- level contributions of type (a) can be computed with standard tree- level methods, 
and compact expressions can be obtained by using MHV-techniques [18]. The one- loop 
terms of type (b) were derived in an analytic form in [19], and form part of the NLO 
corrections to H + 2j final states. The H — > ggg and H — > gqq matrix elements were 
previously known to one loop [20]. In this paper, we compute the two-loop corrections to 
these processes in the heavy top quark limit. We note that expressions for the two- loop 
H —7- ggg helicity amplitudes were previously given in the PhD thesis of one of us. [21] 

The three different contributions must be combined into a parton-level event genera- 
tor program. All three are separately infrared-divergent, and only their sum is finite and 
physically meaningful. To combine the contributions, an infrared subtraction method is re- 
quired. Several methods have been applied successfully in NNLO calculations of exclusive 
observables in the recent past: sector decomposition [22], (^-subtraction [23] and antenna 
subtraction [24-26]. A resulting parton-level event generator will allow an NNLO descrip- 
tion of both H + lj production and of the Higgs boson transverse momentum distribution. 

In addition to their phenomenological importance for precise predictions of collider 
processes, two-loop amplitudes are interesting to investigate fundamental properties of 
quantum field theory at high perturbative orders, aiming to identify regularities, asymptotic 
behaviour and heading towards an all-order understanding of field theory amplitudes. In 
this context, the massless 2—^2 QCD scattering amplitudes at two loops [27] were used to 
determine the high energy Regge trajectories of quarks and gluons at the two-loop level [28]. 
Similarly, the two-loop decay matrix elements for 7* — > qqg [29] and H — > 3 partons fix the 
the two-loop splitting amplitudes [30], which describe the collinear factorization of loop 
amplitudes. 

This paper is structured as follows: in Section 2, we consider the effective coupling of 
the Higgs boson to light partons and define the kinematics. In Section 3, we describe the 
method used to construct the tensor coefficients in the general amplitudes. In Section 4, we 
construct the helicity amplitudes. The derivation of the two-loop corrections to the helicity 
amplitudes, their renormalization and infrared properties are described in Section 5. We 
conclude with an outlook in Section 6. The one-loop and two-loop helicity amplitudes can 
be expressed in compact analytic form, and are enclosed in the Appendix. 

2. Notation and kinematics 
2.1 The effective Lagrangian 

At tree level, the Higgs boson does not couple either to the gluon or to massless quarks. In 
higher orders in perturbation theory, heavy quark loops introduce a coupling between the 
Higgs boson and gluons. As we mentioned in Sec. 1, in the heavy top quark limit, Mt — > 00, 
the H gg coupling becomes independent of Mt. We can therefore integrate out the top quark 
field and formulate an effective Lagrangian, £ e // [7] that couples the scalar Higgs field and 
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the gluon field strength tensor, thereby approximating the Hgg coupling. This large top 
quark mass approximation has been shown to work very well under the condition that the 
kinematic scales involved are smaller than twice the top quark mass [31]. 
The effective Lagrangian reads, 

Ant = — —HG%"G atfiv ■ (2.1) 

where G% u is the field strength tensor of the gluon. The coupling A has inverse mass 
dimension. It can be computed by matching [32,33] the effective theory to the full standard 
model cross section [6]. 

2.2 Kinematics 

We consider the decay of the Higgs boson to three gluons, 

H(p 4 ) — ► gi ( Pl ) + g 2 (p 2 ) +53(P3) , (2.1) 
or into a quark-antiquark pair and a gluon, 

H(p 4 ) — > q( Pl ) + q(p 2 ) + g (P3 ) . (2.2) 
It is convenient to define the invariants, 

S12 = (pi +P2) 2 , S13 = (pi +P3) 2 , s 2 3 = (P2 +P3) 2 , (2.3) 

which fulfill 

P\ = S\2 + S13 + S 2 3 = S123 = Mjj , (2.4) 

as well as the dimensionless invariants, 

X = Si 2 /si23, y = Sl 3 /s 12 3, Z = S 23 /s 12 3, (2.5) 

which satisfy x + y + z = 1. 

3. The general tensors 

The amplitudes \M) can be written as, 

\Mggg) = Sp Up (gi ; g 2 ; g 3 ) £i e 2 el , 

\M qqg ) = Tp(q,q;gy, (3.1) 
while the partonic currents may be perturbatively decomposed as, 

Wsi; 92-,9s) = \V^sf aia2az [s^ p (gi; 92193) + (g) S$ p ( 9l ;g 2 ;93) 

+ {^) 2 S%( gi -g 2 -g3) + 0(al)\ (3.2) 
T p (q; q; g) = Xy/to^T* [t® (q; q; g) + (g) T« (q; q; g) 

+ {^) 2 Tp{q-q-g) + 0{al)], (3.3) 
(i) (i) 

where a s is the QCD coupling constant, and Splp and Tp are the i-loop contributions to 
the amplitude. The SU(3) generators are normalized as tr(T a T b ) = 5 ab /2. 
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3.1 The general tensor for H — > ggg 

The most general tensor structure for the partonic current S /jtl/p (gi; g 2 ; <7 3 ) is given by, 

3 3 

S^u P (gi;g2;g3)4 e 24 = ^2 A ijk p i -eip j -e 2 p k -e 3 + ^2 Bip v eie 2 -e 3 

i,j,k=l i=l 
3 3 

+ ^ CiPi-^ei-^ + ^ D i p i -e 3 e 1 -e 2 

i=l i=l 
= A 211 p 2 -e 1 p 1 -e 2 p 1 -e 3 + A 212 p 2 -e 1 p 1 -e 2 p 2 -e 3 + A 231 p 2 ■ a x p 3 ■ e 2 p\ ■ e 3 

+ A 232 p 2 -t\p 3 -t 2 p 2 -t 3 + A 31 xP3-txpx-t 2 px-t 3 + A 3 i 2 p3-eiPi-e2P2-e3 

+ ^33iP3-eiP3-e2Pi-e 3 + -4 33 2P 3 -eip 3 -e 2 p2-e 3 

+ B 2 e 2 -e 3 p 2 -e 1 + B 3 e 2 -e 3 p 3 -e 1 

+ Ci e 1 -e 3 p 1 -e 2 + C 3 e 1 -e 3 p 3 -e 2 

+ Z?i e 1 -e 2 p 1 -e 3 + £> 2 er£2P2-e 3 , (3.4) 

where the constraints p\ • ei = 0, p 2 • e 2 = and p 3 • e 3 = have been applied. The tensor 
must satisfy the QCD Ward identity when the gluon polarization vectors e\, e 2 and e 3 are 
replaced with the respective gluon momentum, 

(ei ->pi) S lxvp {g 1 ;g 2 ;g 3 )p^e 2 e p 2i = 0, 
(e 2 ~^P2) -> S IJiV p{g 1 ;g 2 ;g 3 )e1p 2 e p 2i = 0, 

(£3 ^P3) -> S lxvp {g 1 ;g 2 ;g 3 )e^e 2 p p 2i = 0. (3.5) 

These constraints yield relations amongst the 14 distinct tensor structures and applying 
these identities give the gauge invariant form of the tensor, 

S^p{gi; 32; SsKe^S = ^21lT 2 n + A311T311 + ^232^232 + A 312 T 312 , (3.6) 

where A^ are gauge independent functions and the tensor structures are given by, 
T 23 2 = P2-eip 3 -e 2 p 2 -e 3 -^ e 2 -e 3 p 2 -ei s 23 - 
T 2 \\ = p 2 -eipi-e 2 pi-e 3 --ei-e 2 pi-e 3 s\ 2 - 
T 3 u = p 3 -e 1 pi-e 2 pi-e 3 -^ei-e 3 pi-e 2 s i3 - 

r 3 i2 = p 3 -e 1 p 1 -e 2 p 2 -e 3 -p 2 -e 1 p 3 -e 2 p 1 -e 3 +-e 1 -e 3 p 3 -e 2 s 12 + -e 1 -e 2 p 1 -e 3 s 23 

1 1 1 1 /nn . 

— ere 3 p r e 2 s 23 + - e 2 -e 3 p 2 -ei s i3 -- e r e 2 p2-e 3 s i3 - - e 2 -e 3 p 3 -ei s± 2 . (3.7) 

The coefficients are functions of the invariants s\ 2 , s 23 and s\ 3 and are further related 
by symmetry under the interchange of the three gluons, 
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The coefficients A^ may be easily extracted from a Feynman diagram calculation 
using projectors such that, 



^2 V(A ijk ) Sf J ,„p(g 1 ;g 2 ;g3)e l teZe% = A ijk , 



(3.9) 



spins 



where the four projectors are given by, 
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Each of the tensor coefficients A; L j k has a perturbative expansion of the form, 



A ijk = X^/Aira 
while the tree-level values are, 



(3.10) 



(3-11) 
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(3.12) 



3.2 The general tensor for H — > qqg 

The most general tensor structure for the partonic current T p (q;q;g) is given by, 

T P (q;q;g)e% = A 1 u(pi)p 3 v(p 2 )pi ■ e 3 + A 2 u{pi)p^v{p 2 )p 2 ■ e 3 + A 3 u(pi)^v{p 2 ), (3.13) 
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where p% ■ €3 = has been applied. The QCD Ward identity yields, 



-43 = ~Pi ■ P3A1 - P2 ■ P3A2 , 
such that the amplitude can be written as, 

T p (q; q; g)el = A 1 (u{p 1 )p 3 v(p 2 )p2 ■ £3 - u{p 1 )^v{p 2 )p2 ■ Pz) 

+A 2 (u{pi)f ? v{p2)Pi ■ £3 - u(Pi)^3 v (P2)Pi ■ P3) (3-14) 
= A 1 T 1 +A 2 T 2 . (3.15) 

The coefficients A{ can be extracted from a Feynman diagram calculation by using 
projectors such that, 

Y^V(A t )T p (q;q;g)e p 3 =A i , (3.16) 

spins 

where the projectors are given by, 

V(M ) ~ (D ~ 2) Tt - — " 4) Tt (3 17) 

v(A\- {D ~ A) tU {d ~ 2) ^ mtt 



Each of the coefficients Ai has a perturbative expansion of the form, 

2 



Ai = A^4 



7ra. s 



(3.19) 



while the tree-level values are simply, 



A f) = 4°) = J_ . (3.20) 

4. Helicity amplitudes 

The general form of the renormalized helicity amplitude \Mggg 2Xa ) for the decay, H{p^) — > 
9i{pi, Ai) + g2(p2, A2) + 53(^3, A3) can be written as, 

IK^ 2A3 > = S w te;92;33)< Al (pi)^ 2 (K)<A 3 (P3) , (4-1) 

where the Aj = ± denote helicity. Similarly, the amplitude for the decay \M. X ^g 2X3 } for the 
decay, H(p 4 ) — > q(pi, Ai) + q(p2, A 2 ) + g(ps,^3) can be written as, 

\M^) = T p (q^;q^;gy 3M (p 3 ) . (4.2) 

The helicity amplitudes can be obtained from the general L>-dimensional tensors of 
Eqs. (3.4) and (3.13) by setting the dimensionality of the Lorentz matrices to be four and 
using standard four-dimensional helicity techniques [34-36]. This corresponds to working 
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in the 't Hooft-Veltman scheme. We use the standard convention of denoting the two 
helicity states of a four-dimensional light-like spinor ^{p) by, 

^±(P) = ^(1±75MP), (4-3) 

with the further notation, 

|p±) = V±(p), (p±|=^±(p). (4.4) 

Particles may thus be crossed to the initial state by reversing the sign of the helicity The 
basic quantity is the spinor product, 

(pq) = {p-\q+), [pq] = [p+\q-], (4.5) 

such that 

(pq)[qp\ = 2p ■ q. (4.6) 

The polarization vector of a outgoing light-like particle with momentum p can then be 
written as 

e^ P ;g) = J q ^ PT) (4.7) 

where q is a light-like reference momentum that satisfies q ■ p ^ but which otherwise can 
be chosen freely. 

Important identities relating spinorial objects are the Fierz rearrangement, 

(p+\l^q+){r + ms+) = 2[pr](sq) (4.8) 

and charge conjugation, 

{p+\r\q+) = {q-\r\p-)- (4.9) 

Substituting Eq. (3.7) into Eq. (3.6), we can express the helicity amplitudes for H — > 
ggg directly in terms of spinor products. It turns out that the only two independent helicity 
amplitudes are \M^ g + ) and |.M++ _ ). The other helicity amplitudes are obtained from 
\M^ + ) and \M^~) by the usual parity relation and by exploiting the symmetry of the 
gluons. Explicitly, choosing pi + \ as reference momentum for we find, 

\M +++ ) = a- ^ , 

999 y/2 (P1 P 2 ) (P2P3 > (PsPl > ' 

\M ++ -) = B- [PlP2]3 (4 10) 

where the coefficients a and (3 are written in terms of the tensor coefficients, 

■S12S13S23 (Si2 S 2 3 , S 13 , „ 

a = oa j4 ^.211 H ^232 ^311 - ^312 

2Mfj \S 2 3 Si 3 S 2 3 

P = S ~fA 211 . (4.11) 
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Likewise (3.14) yields the helicity amplitudes for H — > qqg in terms of spinor products. 
There is only one independent helicity amplitude \-A4q^ + ) and all other amplitudes can 
be obtained from \M q ^ + ) using the usual parity and charge conjugation relations. By 
choosing p\ as reference momentum for 63^3 , we obtain, 



7- 



1 



s/2 \p1P2} ' 

The helicity coefficient 7 is obtained from the tensor coefficients as, 

7 = S12A1 . 



(4.12) 



(4.13) 



As with the tensor coefficients, the helicity amplitude coefficients a, f3 and 7 are vectors 
in colour space and have perturbative expansions, 



(4.14) 



for Q = a, p, 7. The colour factor is T OL = T f} = f a ^ a s and T 7 = If 3 , 



UJ2' 



5. Calculation of the two-loop helicity coefficients 
5.1 Calculation of two-loop Feynman amplitudes 

The calculation of the two-loop Feynman amplitudes contributing to H — > ggg and H — > 
qqg follows closely the calculation of the two-loop helicity amplitudes for 7* — > qqg [29], 
which contribute to the NNLO corrections to e + e _ — > 3j and related event shapes [37,38], 
and of the two-loop helicity amplitudes for qq [39]. We performed two completely 

independent calculations of the amplitudes, which provide a strong internal cross-check on 
the results. 

The Feynman diagrams contributing to the i-loop amplitude \M^) (i = 0, 1,2) were 
all generated using QGRAF [40]. For H — > ggg, there are four diagrams at tree- level, 60 
diagrams at one loop and 1306 diagrams at two loops, while for H — > qqg, we have one 
diagram at tree-level, 15 diagrams at one loop and 228 diagrams at two loops. We use 
dimensional regularization [41-43] with D = 4 — 2e dimensions. We therefore apply the 
.D-dimensional projectors given in Eqs. (3.10) and (3.18) and perform the summation over 
colours and spins using computer algebra methods, mainly implemented in FORM [44]. 
When summing over the polarizations of the external gluons in the projectors, we use the 
axial gauge with a D-dimensional metric. Internal gluons are kept in Feynman gauge, 
resulting in internal ghost contributions to the loop amplitudes. The integrals appearing 
in the individual two-loop diagrams contain up to seven propagators in the denominator, 
and up to five irreducible scalar products in the numerator (i.e. scalar products which can 
not be expressed as linear combinations of the occurring propagators). 

The reduction of the two-loop integrals to a small set of master integrals using integration- 
by-parts (IBP) [45,46] and Lorentz invariance (LI) [47] identities was performed using the 
Laporta algorithm [48], which is based on a lexicographic ordering of the integrals. We 
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used two independent implementations of the Laporta algorithm: the MAPLE and FORM 
based implementation which was developed in the context of [29] and the recently developed 
C++ code REDUZE [49]. Both implementations are based on auxiliary topologies [29], and 
substantial work was required to automate the translation of the momentum assignments 
in the diagrams generated by QGRAF into the momentum sets of the auxiliary topolo- 
gies. This process has been automated in FORM using an iterated shifting and matching 
algorithm for the momenta. 

The two-loop master integrals relevant for this calculation are two-loop four-point 
functions with one leg off-shell. These functions were all computed in [50] in dimensional 
regularization. The results of [50] take the form of a Laurent series in e, starting at e -4 , 
with coefficients containing one- and two-dimensional harmonic poly logarithms (HPLs [51] 
and 2dHPLs [50]), which are a generalization of Nielsen's polylogarithms [52]. Several 
numerical implementations of HPLs and 2dHPLs are available [53] . 

Inserting the master integrals into the amplitudes and truncating the Laurent series 
to the required order, the unrenormalized one-loop and two-loop helicity coefficients are 
obtained. Their Laurent expansion contains HPLs and 2dHPLs up to weight 4. The 
expressions for the master integrals derived in [50] apply to the kinematical situation of a 
1—7-3 decay, while the H + lj production corresponds to a 2 — > 2 scattering process which 
is obtained from the decay kinematics by crossing. The crossing of the amplitudes requires 
the analytic continuation of the master integrals, which is described in detail in [54]. 

5.2 Ultraviolet renormalization 

Renormalization of ultraviolet divergences is performed in the MS scheme. It is carried out 
by replacing the bare coupling «o with the renormalized coupling a s = a s [^i 2 ), evaluated 
at the renormalization scale fi 2 , 

?2 



where 



(5.1) 



S e = (47r) e e- £7 with Euler constant 7 = 0.5772 . . . 

and //q is the mass parameter introduced in dimensional regularization [41-43] to maintain 
a dimensionless coupling in the bare QCD Lagrangian density; fio and p\ are the first two 
coefficients of the QCD /3-function, 

a 11C A -4T R N F 17C 2 A -10C A T R N F -6C F T R N F 

P0 = g , Pi = g , [0.2) 

with the QCD colour factors, 

iV 2 - 1 1 
C A = N, C F = ^ r , T R = -. (5.3) 

The renormalization relation for the effective coupling A is given in [55] as, 



A^ = A 



2 / a -N2 



) +0{o? s ) 



(5.4) 
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We denote the i-loop contribution to the unrenormalized coefficients by n^' U , using 
the same normalization as for the decomposition of the renormalized amplitude (4.14). The 
renormalized coefficients are then obtained as, 



= 5 -i n (i),u _ £^^(o),u 

2e 



3/3o 
2e 

n (2) = 5 -2^(2),U _ ^ S -l n (l),U - ^ - 0(°)' U . 



(5.5) 



For the remainder of this paper we will set the renormalization scale y? = Mfj = S123. 
The full scale dependence of the helicity coefficients are given by, 



Q = AV47ra.(/iWn(°> + fe^ 



£_ 

2ru-- — 1 
V s 123 



SlW + lfoQW In 



+ 



2vr 



2^ 



"^ {2) + ( + ktl (0)> ) In f ^ ) + ^ (0) In 2 f ^ 

V 2 2 / V s 123/ 8 V s 123 



+0(af)|. (5.6) 
5.3 Infrared factorization 

The amplitudes contain infrared singularities that will be analytically canceled by those 
occurring in radiative processes of the same order (ultraviolet divergences are removed by 
renormalization). Catani [56] has shown how to organize the infrared pole structure of the 
one- and two-loop contributions renormalized in the MS scheme in terms of the tree and 
renormalized one- loop amplitudes. This formula for the pole structure is proven [57] from 
the structure of soft and collinear radiation in perturbation theory and can be generalized 
to higher loop order. 

The same factorization of pole terms applies to the helicity coefficients. In particular, 
the infrared behaviour of the one-loop coefficients is given by, 

fiM = /W( e ){i(°) + nW' finite , (5.7) 

while the two-loop singularity structure is, 

+ K) lg)(2e) + flg>( e) ) n<» 

+j£j ) (e)fi (1) + n (2)Jinite , (5.8) 
where the constant K is, 



- 11 - 



For each of the processes under consideration, there is only one colour structure present 
at tree level. Adding higher loops does not introduce additional colour structures, and 
the amplitudes are therefore vectors in a one-dimensional space. Similarly, the infrared 
singularity operators 1$ (e) are lxl matrices in the colour space and are given by, 



on 



2r(i - 



2r(i - e ) 



N i tI + It: ) (Si2 + Si 3 + s 23 ) 



Ne 







4e 2Ne 



(Sl3 + S 2 3)-^(^ + |]S, 



(5.10) 
(5.11) 
,(5.12) 



where, since we have set ji 2 = S123, 



S123 



337 



(5.13) 



Note that on expanding Sjj, imaginary parts are generated, the sign of which is fixed by 
the small imaginary part +i0 of Sij. The origin of the various terms in Eqs. (5.10)-(5.12) 
is straightforward. Each parton pair ij in the event forms a radiating antenna of scale Sij. 
Terms proportional to Sij are canceled by real radiation emitted from leg i and absorbed 
by leg j. The soft singularities 0(1/ e 2 ) are independent of the identity of the participating 
partons and are universal. However, the collinear singularities depend on the identities of 
the participating partons. For each quark we find a contribution of 3/(4e) and for each 
gluon we find a contribution of /3o / (2e) coming from the integral over the collinear splitting 
function. 

Finally, the last term of Eq. (5.8) that involves H^ 2 \e) produces only a single pole in 
e and is given by, 

(5.14) 



4er(l-e) n ' 



(2) 

where the constant H n is renormalization scheme dependent. As with the single pole 
parts of J^(e), the process-dependent can be constructed by counting the number of 
radiating partons present in the event. In our case, 



= H 



(2) 



3H( 2 \ 



(5.15) 



where, in the MS scheme, 



11 q 



9 



7 A 409 

-C3 H 

4 Sd 864 



+ 



96 

vr2 _ 25 \ (N 2 - 1)N F 
48 ~ 216 ) N 



7T 

96 



1 , 5 llvH 

-C3 H 1 

2 Sd 12 144 



N z + 



27 



TT 

72 



89 
108 



, _3 _ 3_ vH 
f 1 2 C3 32 + 8 J N 2 



NN F 



N F 

In' 



(5.16) 
(5.17) 
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At leading order, one can insert the values of the tensorial coefficients given in Eqs. (3.12) 
and (3.20), into Eqs. (4.11) and (4.13) respectively to find, 



a (o) =/3 (o )=7 (o ) = 1 _ (518) 

The renormalized NLO helicity amplitude coefficients can be straightforwardly obtained 
to all orders in e from the helicity coefficients f^ 1 ). For practical purposes they are needed 
through to 0(e 2 ) in evaluating the one-loop self-interference and the infrared divergent 
one-loop contribution to the two-loop amplitude, while only the finite piece is needed for 
the one-loop self-interference. They can be decomposed according to their colour structure 
as follows, 

n {l)Jmite = (n A$> + 1 Bjp + N F C { ^ . (5.19) 

The finite two-loop remainder is obtained by subtracting the predicted infrared struc- 
ture (expanded through to 0(e )) from the renormalized helicity coefficient. We further 
decompose the finite remainder according to the colour casimirs as follows, 



n ^nue = , N , A m + N « B m + i C ( 2 ) + Nf d w + NNFE g) + NlF w 



N 2 n N 



n 



(5.20) 



All one- and two-loop coefficients are given in Appendix A and B respectively. 

To calculate the two-loop contributions to Higgs-boson-plus-jet production at hadron 
colliders, the helicity amplitudes must be crossed to the appropriate kinematical situations. 
Two types crossings are required: 

g(Pi) + 9{P2) -> H(p 4 ) + g(-p 3 ) , q(pi) + q{p2) -»• H(p 4 ) + g(-p 3 ) , (5.21) 
9{P2) + g(ps) -»• H( P4 ) + g(- Pl ) , q(p 2 ) + g(p 3 ) -»• H(p 4 ) + q(- Pl ) . (5.22) 

The definitions of the helicity amplitudes in terms of momentum spinors (4.10,4.12) remain 
unchanged by the crossing, such that only the helicity coefficients a, ft, 7 are to be continued 
to the appropriate kinematical region. The analytical continuation of the polylogarithmic 
functions appearing in two-loop amplitudes is described in detail in [39,54]. We provide the 
one-loop and two-loop coefficients in all relevant analytic continuations in FORM format 
with the arXiv-submission of this article. 



6. Conclusions 

In this paper, we derived the two-loop corrections to the helicity amplitudes for the pro- 
cesses H — > ggg and H — > qqg. Our calculation was performed in dimensional regulariza- 
tion by applying D-dimensional projection operators to the most general tensor structure 
of the amplitude. Our results are expressed in terms of dimensionless helicity coefficients, 
which multiply the basic tree-level amplitudes, expressed in four-dimensional spinors. By 
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applying Catani's infrared factorization formula, we extract the finite parts of the helicity 
coefficients, which are independent on the precise scheme used to define the helicity am- 
plitudes. We provide compact analytic expressions for the two-loop helicity coefficients in 
terms of HPLs and 2dHPLs. 

By crossing the Higgs boson to the final state and two partons to the initial state, 
these amplitudes describe the two-loop corrections to the parton-level process for H + lj 
production and for the transverse momentum distribution of the Higgs boson. To compute 
the NNLO corrections to these processes, the newly derived two-loop amplitudes need to be 
combined with the previously known [18,19] crossings of the tree-level amplitudes for H — > 
5 partons and one-loop amplitudes for H — > 4 partons into a parton-level event generator. 
The tree-level double-real radiation and one-loop real-virtual contributions both contain 
infrared singularities from soft or collinear real radiation. In order to numerically implement 
these contributions, the singular contributions must be subtracted and combined with 
infrared singularities from the two-loop integrals. Up to now, infrared subtraction at NNLO 
has not been fully accomplished for hadron collider processes involving final state jets. In 
the context of dijet production, a first proof-of-principle implementation of the double- 
real radiation contribution exists [26], and the infrared structure of the NNLO subtraction 
terms for hadron collider processes is largely understood [25]. With the amplitudes derived 
in this paper, it should thus become feasible to compute the NNLO corrections to H — > lj 
production and the Higgs boson transverse momentum distribution at hadron colliders. 
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A. One-loop helicity coefficients 

The finite contributions to the renormalized one-loop helicity coefficients, decomposed in 
colour factors according to (5.19) are: 




H(0, 1, z) - H(0, z)G(l -z,y) + G(0, y)H(l, z) - G(0, y)H(0, z)) 
+ G(-z, l-z,y) + G(l, 0, y) - H(l, z)G(-z, y) 

+ H ( _ G (l - z, y) + H(l, z) - H(0, z) - G(0, y)) - ^ 




(A.l) 
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4 1} = o, 



A 



(i) 



I? 



(i) 



c 



(1) 



,4(1) 

7 



S (i) 

7 



C (l) 
^7 



G(l -z,y)- H(l, z) + F(0, z) + G(0, y) 



I(( y + z)(l-y)-z 2 ) 



.7T 

''2' 



(A.2) 



(A.3) 



- G(l - z, 0, y) - (1, 0, z) - G(0, 1 - z, y) - fT(0, 1, z) 

- (0, z)G(l - 2, y) + G(0, y)F(l, z) - G(0, y)tf(0, z)) + G(-z, l-z,y) 
+ G(l,Q,y)-H(l,z)G(-z,y) 

+ ^ ( - G(l - z, y) + (1, z) - ff(0, z) - G(0, y)) 

Z A 1 Z \ 7T 2 .117T 

~ 6(1 - y - z) V 1 - y - z + 1 - y - z J ~ 12 + *~ ' ( j 

, (A.5) 



P ~ 6 



G(l - z, y) - z) + H(0, z) + G(0, y) 



1 - y - z 



1 + 



1 



1 — y — z 1 — y — z 



7T 

>2' 



(A.6) 



3( - G(0, y)tf (1, z) + #(0, z)G(l - z,y) + H(0, 1, z) + G(0, 1 - z, y) 
- G(l, 0, y) + G(l - z, 0, y)) + 6 (tf (1, z)G(-z, y) - G(-z, 1 - z, y)) 
+ 5 (G(0, y) + (0, z)) + y ( - iT(l, z) + G(l :.//)) y ' 



lbr 



— z- 

1 

6 
1 

6 



3( - G(0, y)tf (0, z) + G(l, 0, y) - (1, 0, z)) 
I(-G(0,y)-iT(0,z)) + 2(^(1, z)-G(l-z,y)) + 



2y 



vr 2 27 3(z - 1) 
_ _ _ 

10 



2y 

.vr 



(A.7) 
,(A.8) 
(A.9) 



B. Two- loop helicity coefficients 

The finite contributions to the renormalized two-loop helicity coefficients, decomposed in 
colour factors according to (5.20) are: 



G {1 _ Zj 1-^,0, 3/)- G(l - z, -z, 0, 1 - z, y) + G(l - z, 1 - z, 0, 0, y) 



+ G(l - z, 0, -z, 1 - z, y) + G(l - z, 0, 1 - z, 0, y) - G(l - z, 0, 1, 0, y) 
+ G(l - z, 0, 0, 1 - z, y) + (1, 1,0, 0, z) + 0, 1, 0, z) + H(l, 0, 0, 1, z) 
+ H(l, 0, 0, z)G(l - z, y) - F(l, 0, z)G(-z, 0, y) + #(1,0, z)G(l - z, -z, y) 
- G(l, 0, y)H(l, 0, z) + H(l, z)G(l - z, -z, 0, y) + z)G(l - z, 0, -z, y) 
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- H(l, z)G(l - z, 0, 0, y) - G(0, -z, l-z,y) + G(0, -z, y)H(l, 0, z) 

- G(0, l-z,-z,l- z, y) + G(0, 1 - z, 1 - z, 0, y) + G(0, 1 - z, 0, 1 - z, y) 
+ G(0, 1 - z, -z, y)tf (1, z) - G(0, -z, 1 - z, y) + G(0, 1 - z, 1 - z, y)H(0, z) 

- G(0, 1 - z, 1, 0, y) - G(0, 1 - z, 0, y)H(l, z) - G(0, 1 - z, 0, y)tf (0, z) 

- fT(0, z)G(l - z, 0, 1 - z, y) - H(0, z)G(l - z, 0, 0, y) - G(0, y)tf(l, 1,0, z) 

- G(0, y)tf(l, 0, 1, z) + G(0, y)tf (1, 0, 0, z) - G(0, 1, 1 - z, 0, y) + (0, 1, 1, 0, z) 

- H(0, 1, 1, z)G(0, y) - G(0, 1, 0, 1 - z, y) + H(0, 1, 0, 1, z) 

+ #(0, 1, 0, z)G(l - z, y) + F(0, 1, 0, z)G(0, y) + G(0, 1, 0, y)H(l, z) 

- G(0, 1, 0, y)H(0, z) - H(0, 1, z)G(-z, 0, y) - H(0, 1, z)G(l, 0, y) 

+ fl-(0, 1, z)G(0, -z, y) + G(0, 0, 1 - z, 1 - z, y) - G(0, 0, 1 - z, z) 

- G(0, 0, 1 - z, y)iJ(0, z) + F(0, 0, 1, 1, z) + fT(0, 0, 1, z)G(0, y) 

+ il(0, 0, z)G(l - z, 1 - z, y) + i7(0, 0, z)G(l - z, 0, y) + i7(0, 0, z)G(0, l-z,y) 

+ G(0, 0, y)i7(l, 1, z) - G(0, 0, y)H(l, 0, z) - G(0, 0, y)H(0, 1, z) 

+ G(0, 0, y)F(0, 0, z) + fT(0, z)G(-z, 1 - z, 0, y) + F(0, z)G(-z, 0, 1 - z, y) 

- il(0, z)G(l - z, -z, 1 - z, y) - i?(0, z)G(l - z, 1 - z, 0, y) 

- H(0, z)G(l - z, 1, 0, y) + H(0, z)G(l, 1 - z, 0, y) + H(0, z)G(l, 0, 1 - z, y)) 

+ ( - G(-z, 1 - z, 1 - z, 0, y) - G(-z, 1 - z, 0, 1 - z, y) 

- G(-z, 0, 1 - z, 1 - z, y) - G(l - z, 1 - z, -z, 1 - z, y) + G(l - z, 1 - z, 1, 0, y) 
+ G(l - z, 1, 1 - z, 0, y) + G(l - z, 1, 0, 1 - z, y) - G(l - z, 1, 0, 0, y) 

- G(l - z, 0, -z, 1 - z, y) - G(l, 1 - z, 0, 0, y) - 1, 0, z)G(-z, y) 

- 1, z)G(-z, 0, y) - G(l, 0, 1 - z, 0, y) + 0, 1, z)G(-z, y) 

- G(l, 0, 0, 1 - z, y) + G(l, 0, 0, y)tf(l, z) - H(1,0, z) + (1, 0, z)G(-z, 1 - z, y) 

- 0, z)G(l - z, 1 - z, y) + H(l, 0, z)G(l - z, 0, y) + H(l, z)G(-z, 1 - z, 0, y) 
+ #(1, z)G(-z, 0, 1 - z, y) + z)G(l - z, 1 - z, -z, y) 

- z)G(l - z, 1, 0, y) - G(0, -z, 1 - z, 1 - z, y) + G(0, -z, 1 - z, 0, y) 
+ G(0, -z, 0, 1 - z, y) + G(0, -z, 1 - z, y) - G(0, -z, y)il(l, 1, z) 

+ G(0, -z, 1 - z, y)tf(0, z) - G(0, -z, 0, y)F(l, z) + G(0, -z, 1 - z, y)F(l, z) 
+ H(0, 1, 1, z)G(-z, y) - H(0, 1, z)G(-z, 1 - z, y) + H(0, 1, z)G(l - z, 1 - z, y) 
+ H(0, 1, z)G(l - z, 0, y) + G(0, 0, -z, 1 - z, y) - G(0, 0, -z, y)F(l, z) 

- G(0, 0, 1, 0, y) - tf(0, z)G(-z, 1 - z, 1 - z, y) + H(0, z)G(l - z, 0, 0, y) 
-iJ(0,z)G(l,0,0,y)) 

+^ ( - G(0, -z, 1 - z, y)iJ(0, z) - i?(0, 1, z)G(l - z, -z, y) 
-iJ(0,0,l,z)G(l-z,y)) 

+2( - G(-z,-z,-z, 1 - z,y) + G(-z, -z, 1 - z, 1 - z,y) 
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+ G(-z, 1 - z, -z, 1 - z, y) + G(l - z, -z, -z, l-z,y)- G(l, 1, 1, 0, y) 

+ G(l, 1, 0, 0, y) + H(l, 1, 0, z)G(l -z,y) + H(l, 1, z)G(-z, -z, y) + G(l, 0, 1, 0, y) 

+ (1, z)G(-z, -z, -z, y) - H(l, z)G(-z, -z, 1 - z, y) 

- H(l, z)G(-z, 1 - z, -z, y) - H(l, z)G(l - z, -z, -z, y) + G(0, 1, 1,0, y) 
+ H(0, 1, z)G(-z, -z, y) + H(0, 0, 1, z)G(-z, y)) 

+Y A { - G (°> 1 - ^ (°. *) + H (°> !. ^ G (°. y) - # (°> ^M 1 - °. y) 

+ G(0,y)tf(l,0,z)) 

+ ¥ ( " G( ~ z ' ~ z ' 1 " z ' y) + G(1 " z ' 1 " z ' °' y) + G(1 " z ' °' 1 " z ' y) 

+ G(l - z, 0, 0, y) - G(l, 1, 0, y) + H(l, 0, 0, z) - H(l, 0, z)G(l - z, y) 
+ H(l, z)G(-z, -z, y) - H(l, z)G(l - z, 0, y) + G(0, l-z,l-z,y) 

- G(0, 1 - z, y)H(l, z) + G(0, 1 - z, 0, y) - fT(0, 1, 1, z) + fT(0, 1, z)G(-z, y) 

+ G(0, 0, 1 - z, y) + #(0, 0, z)G(l - z, y) + #(0, 0, z)G(0, y) - G(0, 0, y)H(l, z) 
+ G(0, 0, y)H(0, z) + fT(0, z)G(l -z,l-z,y) + G(0, y)F(l, 1, z)) 

+^ ( - H(l, 0, 1, z) + F(0, 1, z)G(l -z,y) + H(0, 1, 0, z)) 

+y ( - G(l, 0, 0, y) - G(-z, 1 - z, 1 - z, y) - H(l, 1, z)G(-z, y) 

+ H(l,z)G(-z,l-z,yj) 

55 / \ 

2 

+y (G(l, 1, y) - G(0, 1, y) + (1, 1, z) - G(l, 0, y) + (0, 0, z) + G(0, 0, y)) 

+ (f _ nr) ( " G(1 " z ' y) + z) " ^ (0 ' z) " G(0 ' y) ) 

+ ( " if + y ) ( " G(1 " Z; 1 " z ' v) ~ H (1, lj z) + z)G(1 " *' y) 

-#(0,0,z)-G(0,0,y)) 

/49 7T 2 \ / 

+ - -g- ) ( - G(l - z, 0, y) - G(0, 1 - z, y) - (0, z)G(l - z, y) 
+ G(0,y)tf(l,z)-G(0,y)tf(0,z)) 

+ (l§ + ( G( " Z ' 1 " Z ' V) + G(1> °' y) " Z)G( " Z ' y) ) 
/245 vr 2 \ , , /389 vr 2 \ TW 

/13 451tt 2 \^. . /13 I2657r 2 \ rr/1 , 

- (t + ^ir) G<1 - ~- v) + ( s + ^r) H ^ z) 

/13 1133tt 2 \/ rr/ . ^. x \ llvr 2 G(l,y) 
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1 /5029vr 2 ^ : 99C 3 , 3vr 4 1321^ 



36 V 24 



- 72C 4 + 



4 + 16 



+ 1 ((y + z)(l - y) - z 2 ) + G(l, 0, y) + G(-z, l-z,y)- H(l, z)G(-z, y) 
+ G(0,y)H(l,z)-G(0,y)H(0,z) 

+1 ( - G(l - z, 0, y) - H(l, 0, z) - G(0, l-z,y)- H(0, 1, z) 
- if (0, z)G(l -z,y)- G(0, y)H(l, z) + G(0, y)H(0, zj) 

-li( G (i-,,,)-^(M)) + ^(o,,) 3 ^ 2471 

/25z 1 15z 2 



2 + 18 



f ll2 ( l+ l-y-z ) 4(1 -y-z) 



2z 



yz 5z 2 . 
6 6 1 — y — z 



+ 



(1 - y- z) 1 

/ 25y _ 9y yz_ _ lhy 2 
V12z 4 6 \z 



G(0,y)H(0,z)+H(l,0,z)-G(l,0,y) 
+ ^ 2 + ^ + >-2 y + y 2 )))x 



G(l - z, 0, y) - G(l, 0, y) - G(-z, l-z,y) + H(l, z)G(-z, y) + G(0, 1 - z, y) 

+ tf(0,l,z)-G(0,y)tf(l,z) 

/25z 9z yz lSz^ 9 5z^,2z I 9 ,,\ 
+ (Wy-^ + \-^ + Z + ^ + ^ {l - 2Z + Z)) ) X 

H(0, z)G(l - z, y) - G(-z, 1 - z, y) + H(l, z)G(-z, y) 

+ ^(63 - 93(y + z)+4yz + — (l-2z + z 2 ) + ^(l-2y + y 2 ) 



+ 30(y 2 + z 2 )) 



G(l-z,y)-H(l,z) 



4 f - 63z + 60z 2 - 30z 2 (l - z)( - + - ) 

36 V y \ — y — z 

l/93z(l-z) 145y 27yz 79y 2 30z 

36 V 2 2~ + 2 + ~T~ + 
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7T 
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7z 2 19yz 172/(1 - y) 
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5z 2 



1 - y - z 



9 



+ 26y(l-y-z))tf(0,z) 



( - 1 + 2z - z 2 ) 



2 V36(l -y-z) 9 ' 18(1 - y - z) ^ Z 2 + l-y-z 



(1 -2z + z 2 )) 



36 



36 
55 / 



+ 



36 



— [ - H(0, z)G(l - z, y) - H(0, z)G(0, y) - (0, 1, z) 
- 2#(1, z)G(-z, y) + (1, z)G(0, y) - H(l, 0, z) - G(l - z, 0, y) 
+ 2G(-z, 1 - z, y) - G(0, 1 - z, y) + 2G(1, 0, y) + H( _ # ( 0) z ) + z) 
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(B.l) 



4 2) 



77vr 2 t 3( 3 , 185" 

, (B.2) 

, (B.3) 

12 U 2 ' - - ■ z 

+ G(0, 1 - z, y) + F(0, 1, z) - G(0, y)# (1, z) - G(l, 0, y) + z)G(-z, y) 



■ ( J,( - 1 + 2y - y 2 ) + -{I - y 2 ) - 4 + 2z + y) 



+ G(i-z,o,y) 



- G(-z,l-z,y) 

z ( 
+ T2( 



( _ i + 2z - z 2 ) + 

(l-y-z) 2 1-y- 2 



(1 - z 2 ) - 2 - z - 2y) 



fl-(l,0,z) + G(0,y)if(0,z) -G(l,0,y) 

z ( z , _ , 2 



+ 



12 V 



; -( - 1 + 2z - z 2 ) + -(1 - z 2 ) - 4 + z + 2y) 



tf(0,z)G(l-z,y) 



+ H(l,z)G(-z,y)-G(-z,l-z,y) 
1/15 3z 



+ 



36 V 2 y 



+ -(-l + 2z-z 2 ) + ^(-l+2y-y 2 )-?(y 2 + z 2 ) 



-3(y + z + yz)) 



G(l-z,y)-F(l,z) 

l(3z 2 (l-z)(I + T3 l_)-6z-l|!)^0,z) 

3y 2 (l - y) 



36 v y 
1 / 3z 



36 V 1 — y — z 

,2 



(l-2z + z^)-3z(l-z) + 



201 

-r + l8Cs 



+ 



6y-3yz-^|-)G(0,y)-l 

7T 2 / Z Z 2z 

s (— + 



y - z 



z 1 — y — z 



2 2 Z 3 

- 2z 2 - 



7T 

+ *4' 



i((y + z)(l-y)-z 2 ) 

)-z(2 + z + 2y)) 
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1-y- z' 



l - ( + G(0, 1 - z, y)tf (0, z) - (0, 1, z)G(0, y) + iJ(0, z)G(l - z, 0, y) 



-G(0,y)#(l,0,z) 
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+ i ( - G(l - z, 0, y) - ff(0, z)G(l - z, y) + G(0, y)i7(l, z) - G(0, y)H(0, z)) 
103 41vr 2 \^ /i , /103 77r 2 \ rr/1 , /35 5vr 2 \ , 



+ { 18 + 8 



/ 103 _ 5t^_ 



)G(0,y)-^-(G(l,y) 



27G(l-z,y)\ 1 /1781 63C 3 



+ 1 J ~ 6l~L2- + 



! )1 



+ 



1879^ 
24 

+ ^ ((y + *)(i - y) - 6( - G(-z, 1 - z, y) - G(l, 0, y) 
tf(l,z)G(-z,y)) 

+ 3( - G(0, y)H(l, z) + G(l - z, 0, y) + (1, 0, z) + G(0, 1 - z, y) 
+ H(0, 1, z) + if(0, z)G(l - z, y)) + G(0, y)#(0, z) + 2G(0, y)#(0, z) 

7T 2 275" 

+ 23G(1 -z,y)- 23#(1, z) - 19JJ (0, z) + 



-7y. 1, lly 3 z 2 1 . 3y 2 

1 - - + - + - - 2 + 2z + y - -2- 

b z I2z z y y z 



G(l-z,0,y) 



+ G(-z, 1 - z, y) + G(l, 0, y) - z)G(-z, y) - G(0, 1 - z, y) 
-H(Q,l,z) + G(Q,y)H(l,z) 

f-lz. 1, llz 3 y 2 1 . 3z 2 \ 

+ (-«- 1 " - + TTT (- + -"2 + 2y + z 

V 6 y 12y 2 z z y / 

- tf(l,z)G(-z,y) - fT(0,z)G(l - z,y) 

llz 3 

y — z' 12(1 — y — z)^ 



G(-z, 1 - z,y) 



+ 1- 1 ) + 

V b 1 — y — z 

3z 2 \ 

+ 2 (l-y-z) + ,)- T -^) 



i( (l^)! + I_ 2 



^(l,0,z) + G(l,0,y)-G(0,y)^(0,z) 



- 20 - 



1 / 57 33z, 1 n 
36 V 2 y 

7yz 



33y, 



123. 

+ — iv + z) 2 



G(l-z,y)-tf(l,z) 



+ 



77?/ 



If 

36 V 

(-l + y + z))tf(0,z) 



« (1 + tf) + »£ + 3 to »(i-,)(I + T -l_; 



1 / 33z . 2n 105z 2 i 77 27z 33y 



-43y- 



33y 



7T 2 /7z 



. ))G(0,y) 

77. 

(1 



) + 



llz 2 



2 V18 v 1-y-z' 36(l-y-z) 



: (-l-(l-y-z) 2 + 2z 



2z{l-y-z)-z z ) + - 



y - z 



+iir 



^ (H(Q, z)G(1 -z,y) + H(0, z)G(0, y) + H(0, 1, z) 



+ 2iT(l, z)G(-z, y) - H(l, z) * G(0, y) + #(1, 0, z) + G(l - z, 0, y) 

- 2G(-z, l-z,y) + G(0, 1 - z, y) - 2G(1, 0, y) + H ( + H (0, z) - H(l, z) 



+ G(l - z, y) + G(0, y))) + !g (z(z - 1) - y(l - y - z)) + ^ - g 

^ 3(G(1 - z, 1 - z, y) + H(l, 1, z) - (1, z)G(l - z, y) + fT(0, 0, z 
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+ H(l, 0, 0, z)G(l -z,y)- H(l, 0, z)G(-z, 0, y) + H(l, 0, z)G(l - z, -z, y) 

- G(l, 0, y)H(l, 0, z) + H(l, z)G(l - z, -z, 0, y) + H(l, z)G(l - z, 0, -z, y) 
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- G(0, l- z ,-z,l- z, y) + G(0, 1 - z, -z, y)H(l, z) + G(0, 1 - z, 1 - z, y)F(0, z) 
+ G(0, 1 - z, 0, 1 - z, y) - G(0, 1 - z, 1, 0, y) + G(0, 1 - z, 1 - z, 0, y) 

- G(0, 1 - z, 0, y)tf (1, z) - G(0, 1 - z, 0, y)#(0, z) - G(0, 1, 1 - z, 0, y) 

+ fT(0, 1, 1, 0, z) - #(0, 1, 1, z)G(0, y) - G(0, 1, 0, 1 - z, y) + H(0, 1, 0, 1, z) 
+ #(0, 1, 0, z)G(l -z,y) + H(0, 1, 0, z)G(0, y) + G(0, 1, 0, y)H(l, z) 

- G(0, 1, 0, y)iJ(0, z) - fT(0, 1, z)G(-z, 0, y) - H(0, 1, z)G(l, 0, y) 

+ fl"(0, 1, z)G(0, -z, y) + fl-(0, z)G(-z, 1 - z, 0, y) + H(0, z)G(-z, 0, 1 - z, y) 

- fT(0, z)G(l - z, -z, 1 - z, y) - fT(0, z)G(l - z, 1, 0, y) + fT(0, z)G(l, 1 - z, 0, y) 
+ fT(0, z)G(l, 0, 1 - z, y) - fT(0, z)G(l - z, 1 - z, 0, y) 

- fl-(0, z)G(l - z, 0, 1 - z, y) + F(0, z)G(l - z, 0, 0, y) - G(0, y)F(l, 1, 0, z) 

- G(0, y)il(l, 0, 1, z) + G(0, y)iJ(l, 0, 0, z) + G(0, 0, 1 - z, 1 - z, y) 

- G(0, 0, 1 - z, y)il(l, z) - G(0, 0, 1 - z, y)iJ(0, z) + i7(0, 0, 1, 1, z) 

+ fl-(0, 0, 1, z)G(0, y) + fT(0, 0, z)G(l - z, 1 - z, y) + H(0, 0, z)G(l - z, 0, y) 
+ H(0, 0, z)G(0, 1 - z, y) - G(0, 0, y)F(l, 0, z) + G(0, 0, y)i7(l, 1, z) 

- G(0, 0, y)F(0, 1, z) + G(0, 0, y)H(0, 0, z) 

+ 2 ( - G(-z, 1 - z, 1 - z, 0, y) - G(-z, 1 - z, 0, 1 - z, y) 

- G(-z, 0, 1 - z, 1 - z, y) - G(l - z, 1 - z, -z, 1 - z, y) + G(l - z, 1 - z, 1, 0, y) 
+ G(l - z, 1, 1 - z, 0, y) + G(l - z, 1, 0, 1 - z, y) - G(l - z, 1, 0, 0, y) 

- G(l, 1 - z, 0, 0, y) - 1, 0, z)G(-z, y) - 1, z)G(-z, 0, y) 

- G(l, 0, 1 - z, 0, y) + 0, 1, z)G(-z, y) - H(l, 0, z)G(l - z, 1 - z, y) 
+ H(l, 0, z)G(l - z, 0, y) - G(l, 0, 0, 1 - z, y) + G(l, 0, 0, y)H(l, z) 

+ H(l, 0, z)G(-z, 1 - z, y) + z)G(-z, 1 - z, 0, y) + z)G(-z, 0, 1 - z, y) 
+ z)G(l - z, 1 - z, -z, y) - #(1, z)G(l - z, 1,0, y) - G(0, -z, 1 - z, 1 - z, y) 
+ G(0, -z, 1 - z, 0, y) + G(0, -z, 0, 1 - z, y) - G(0, -z, 0, y)F(l, z) 
+ G(0, -z, 1 - z, y)F(l, z) - G(0, -z, y)il(l, 1, z) + i7(0, 1, 1, z)G(-z, y) 

- #(0, 1, z)G(-z, 1 - z, y) + H(0, 1, z)G(l - z, 1 - z, y) + il(0, 1, z)G(l - z, 0, y) 
+ G(0, 0, -z, 1 - z, y) - G(0, 0, -z, y)F(l, z) - H(0, z)G(-z, 1 - z, 1 - z, y) 
- J ff(0,z)G(l,0,0,y) -G(0, 0,1,0, y)) 

+ 3( - #(0, 1, z)G(l - z, -z, y) - fT(0, 0, 1, z)G(l - z, y)) 



+ 4^ - G(-z, -z, -z, 1 - z, y) + G(-z, -z, 1 - z, 1 - z, y) 
+ G(-z, 1 - z, -z, 1 - z, y) + G(l - z, -z, -z, 1 - z, y) - G(l, 1, 1,0, y) 
+ G(l, 1, 0, 0, y) + H(l, 1, 0, z)G(l - z, y) + H(l, 1, z)G(-z, -z, y) + G(l, 0, 1, 0, y) 
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+ H(l, z)G(-z, -z, -z, y) - H(l, z)G(-z, -z, 1 - z, y) 

- H(l, z)G(-z, 1 - z, -z, y) - H(l, z)G(l - z, -z, -z, y) + G(0, 1, 1, 0, y) 
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+ 
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- H(0, z)G(l - z, 0, y)) - G(l - z, 1, 0, y) - G(l, 1 - z, 0, y) - 1, 0, z) 

- G(l, 0, 1 - z, y) + G(l, 0, y)H(l, z) - H(l, 0, z)G(l - z, y) 
+ H(0,l,0,z) + H(0,1, z)G(l - z,y) 

+ 2 ( - G(-z, -z, 1 - z, y) + G(l - z, 1 - z, 0, y) + G(l - z, 0, 1 - z, y) 
+ G(l - z, 0, 0, y) - #(1, 0, 1, z) + #(1, 0, 0, z) + #(1, z)G(-z, -z, y) 
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+ G(0, 0, y)#(0, z) - G(0, 1 - z, y)#(l, z)) 

+ 3( - G(l - z, -z, 1 - z, y) + H(l, z)G(l - z, -z, y)) 
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- G(l, 0, y)#(l, z) - H(0, 1, z)G(l - z, y) 

+ 2 ( + G(-z, -z, i- z , y )-G(l-z,l- z, 0, y) - G(l - z, 0, 1 - z, y) 

- G(l - z, 0, 0, y) - H(l, 1, 0, z) + z)G(l - z, 0, y) + H(l, 0, 1, z) 

- H(l, 0, 0, z) - H(l, z)G(-z, -z, y) - G(0, 1 - z, 1 - z, y) - G(0, 1 - z, 0, y) 
+ G(0, 1 - z, y)H(l, z) - H(0, 1, z)G(-z, y) + (0, 1, 1, z) - G(0, 0, 1 - z, y) 

- H(0, 0, z)G(l - z, y) - (0, 0, z)G(0, y) + G(0, 0, y)H(l, z) - G(0, 0, (0, z) 

- if (0, z)G(l - z ,\-z,y)- G(0, y)#(l, 1, z)) 

+ 3(G(1 - z, -z, 1 - z, y) - H(l, z)G(l - z, -z, y) - if (0, 1, 0, z)) 
+ 4(G(-z, 1 - z, 1 - z, y) + ii(l, 1, z)G(-z, y) + G(l, 0, 0, y) 

- if (1, z)G(-z, 1 - z, y) + G(0, 1, 0, y)) 



1 

36 



20(G(-z, 1 - z, y) - H(l, z)G(-z, y) + G(l, 0, y)) 

- 21 ( - if (0, 1, z) - ff(l, 0, z)) + 33 ( - G(l - z, 1 - z, y) - H(l, 1, z) 

+ if (1, z)G(l - z, y) - if (0, 0, z) - G(0, 0, y)) - G(l - z, 0, y) - G(0, 1 - z, y) 

- if (0, z)G(l - z, y) + G(0, y)ff (1, z) - G(0, y)H(0, z) + - (G(l - z, y) 

- H(l, z)) - 9G(0, y)if (1, 0, z) + -^ff (0, z) + ^G(0, y) 



- 26 - 



+ 



yz 



G(l - z, 0, y) + G(0, 1 - z, y) + H(0, 1, z) 



12(1- y _ z )2 

+ if(0, z)G(l - z, y) - G(0, z) + 2 ( - G(-z, 1 - z, y) + z)G(-z, y) 

2 n 

- r-^ H(l, z)) + ( - G(0, 1, 0, y) - G(0, y)#(l, 0, z)) 

1 — y — z / 1 — y — z V / 



+ -(G(l-z,y)-g(l,z))- 



12ff(l ,l,0,z) _ 7r 2 G(0,y) + 12C 3 

y — z 1 — y — z 1 — y — z 



+ 



11tt 2 z 2 374 
(1-y-z) 2 "T 

z 

+ 



+ 3tt" 



:!Gf , _ n . v 3(l - y - z) 2 + 3z(l -y-z) + 2z 2 ) 6(g(1, 1, 0, y) 
-#(1,1,0,*)) +3(-F(0,l,0,z) -G(0,l,0,y) - (0, z)G(l, 0, y) 
- G(0, y)F(l, 0, z)) + . 2 ( - G(l, y) - z) + ^ - 

+ 36(1-^-^)3 ( 3 (! " y " ^) 2 " 3(1 - y - z) + 2) [ 6 (F(1, 1, 0, z) - Cs) 
+ 3(G(0, 1, 0, y) + G(0, y)il(l, 0, z)) + 7r 2 iJ(l, . 



z) + 



vr 2 G(0,y) 



36z V z 



- G(l - z, 0, y) + G(-z, 1 - z, y) - (1, z)G(-z, y) 
- <-*'((). ! - :-.(/)- //((), l,z) + G(0,y)^(l,z)+G(l,0,y) 

+ G(-z, 1 - z, y) - (1, z)G(-z, y) - (0, z)G(l - z, y) 



z / 33z 

H 42 H 

36y V y 



11 /y z 
+ —(- + - 

12 Vz y 



+ 



G(l-z,y) + H(l,z) 
12 



(- 



36(1 — y — z) 2 V z(l — y — z) z 
186 198z 



_ 15 + 39(1 - y - z) + 



99 



1 - y - z 



+ 141 - 



198z 
+ z h 



z 

99z 2 



(1-y-z) 2 1-y-z {l-y-zf 



(1-y-z) 2 
-#(1,0,*) 



+ G(l, 0, y) - G(0, y)tf (0, z) - 1 V Z - G(0, y) 



( - 6 + 6(1 - y - z) - 3z 2 ) H(l, 0, z) 



3 6(1 _ 9 _ z)2 ( -12-3^(0,^(0,;) 

1 / 1 xx 33, 1 , 1 ,87 

36\(1 — yY z 1 — y — z [\ — y — zY 2 



1 


36(1 


- y - z) 2 




1 


36(1 


- y - z) 2 




z 



(i-y) 

+ 15(1 -y-z) + 



173z 125z(l -y-z) AO _ 2s 33y 



43z 2 )--^)G(0,y) 



-27- 



1/33, 1 



36 V y 1 — z 
1 



39 



+ 



33 



+ 



+ 



(i-y-z) 

99z 3 



(1 - z) 2 ' (1 - z)(l -y- z) ' 1 - z 

, 25z 131z(l- y-z) 99z 2 
2 (33(1 - y - z) + — + ^ » >- 



1-y-z 



+ U2z 2 



1-y-z 

1 /1115 
~ 36 V 12 



H(0,z) 



+ 39(-^— + -^—) + 



7T 



1-y • 1-z' • (l-y-z)*( 1 



45(1 -y-z) 



+ (1-y-z) - 



2z 



2vr 2 

33z 2 31z 2 



13z 9z(l-y-z) 
{l-y-z) + ~2 2 + 2(1 - y- z) 2 1 



+ 12z 2 



33z 4 



2(l-y-z) 2 ' 

^ ( + H(0, z)G{l -z,y) + H(0, z)G(0, y) + H(0, 1, z) 
+ 2H(1, z)G(-z, y) - H(l, z)G(0, y) + H(l, 0, z) + G(l - z, 0, y) 
- 2G(-z, l-z,y) + G(0, l-z,y)- 2G(1, 0, y) + y ( + (0, z) - H(l, z) 



+ G(l-z,y) + G(Q,y)) 



65yz 



+ 



7tt 2 



71 



,(2) 



72(1 -y-z) 2 ' 144 18 
^ (g(1 -z,l-z,y) + H(l, 1, z) - #(1, z)G(l - z, y) + fT(0, 0, z) + G(0, 0, y)) 



(B.11) 



+ ( + G(l - z, 0, y) - 0, z) + G(0, 1 - z, y) - H(0, 1, z) 
+ if(0, z)G(l - z, y) - G(0, y)tf (1, z) + G(0, y)#(0, z)) 
+ ^ ( - G(0, y) - (0, z) + H(l, z) - G(l - z, y)) 



yz 



36(1 - y- z) 2 

.57T 



( - G(0, y) - (0, z) + 3H(1, z) - 3G(1 - z, y) + y 



36 



- H(0, z) + H(l, z) - G(l - z, y) - G(0, y) + 2 + 



10 

3 

yz 



A? 



G(l-z,0,0,y) 1 
3 2 



(1-y-z) 2 J 
( - G(l - z, 1, 1, 0, y) + G(l, 1 - z, -z, 1 - z, y) 



297T 2 

144 
,(B.12) 



- G(l, 1 - z, -z, y)tf (1, z) + G(l, 1 - z, 1, 0, y) + H(l, 1, 0, 0, z) 

- G(l, 0, -z, y)il(l, z) - G(l, 0, 1,0, y) + 0, 1, z)G(l, y) + G(l, 0, -z, 1 - z, y) 

- 0, 0, z)G(l - z, y) + G(l, 0, y)H(l, 0, z) + G(0, -z, 1 - z, y)F(0, z) 

- G(0, -z, y)tf (1, 0, z) + G(0, 1 - z, -z, 1 - z, y) - G(0, 1 - z, -z, y)tf(l, z) 
+ G(0, 1 - z, 0, y)H(0, z) + G(0, 1 - z, y)H(l, 0, z) + G(0, 1, 0, y)tf (0, z) 

- #(0, 1, z)G(l, 1 - z, y) - F(0, 1, z)G(0, -z, y) + H(0, 1, z)G(0, 1 - z, y) 
+ #(0, 0, 1, z)G(l,y) + H(0, 0, 1, z)G(0, y) - H(0, 0, z)G(l - z, 0, y) 

+ G(0, 0, y)tf (1, 0, z) + G(0, 0, y)H(0, 0, z) + (0, z)G(l - z, 1, 0, y) 

- fT(0, z)G(l - z, 0, 0, y) - H(0, z)G(l, 1 - z, 0, y) - F(0, z)G(l, 0, 1 - z, y) 



- 28 - 



+ G(0, y)H(l, 0, 1, z) + G(0, y)tf(l, 0, 0, z)) 

+ | ( - G(l, 1 - z, 0, y) - G(l, 0, 1 - z, y) + G(l, 0, y)tf (1, z) 

-tf(0,0,z)G(0,y)) 

+ G(l - z, -z, 1 - z, 0, y) + G(l - z, -z, 0, 1 - z, y) - G(l - z, 1 - z, 0, 0, y) 

- G(l - z, 1, 1 - z, 0, y) - G(l - z, 1, 0, 1 - z, y) + G(l - z, 0, -z, 1 - z, y) 

- G(l - z, 0, 1 - z, 0, y) + G(l - z, 0, 1, 0, y) - G(l - z, 0, 0, 1 - z, y) 

+ G(l, 1 - z, 0, 0, y) + G(l, 1, 1, 0, y) - G(l, 1, 0, 0, y) + G(l, 0, 1 - z, 0, y) 
+ 0, 1, z)G(l - z, y) + 0, 0, z)G(-z, y) + G(l, 0, 0, 1 - z, y) 

- #(1, 0, 0, z)G(-z, y) - G(l, 0, 0, y)il(l, z) + 0, z)G(-z, 0, y) 

- H(l, 0, z)G(l - z, -z, y) - (1, z)G(l - z, -z, 0, y) + z)G(l - z, 1, 0, y) 

- H(l, z)G(l - z, 0, -z, y) + z)G(l - z, 0, 0, y) - G(0, 1 - z, 1 - z, 0, y) 

- G(0, 1 - z, 1 - z, y)tf (0, z) + G(0, 1 - z, 1, 0, y) - G(0, 1 - z, 0, 1 - z, y) 

+ G(0, 1 - z, 0, y)F(l, z) + iJ(0, 1, 1, 0, z) - G(0, 1, 1,0, y) + H(0, 1, 1, z)G(0, y) 

- H(0, 1, 0, z)G(l - z, y) + F(0, 1, z)G(-z, 0, y) - tf(0, 1, z)G(l - z, 0, y) 

- G(0, 0, 1 - z, 1 - z, y) + G(0, 0, 1 - z, y)F(l, z) - #(0, 0, 1, 1, z) 

- tf(0, 0, z)G(l - z, 1 - z, y) - G(0, 0, y)tf(l, 1, z) - (0, z)G(-z, 1 - z, 0, y) 

- iJ(0, z)G(-z, 0, 1 - z, y) + H(0, z)G(l - z, -z, 1 - z, y) 
+ H(0, z)G(l - z, 1 - z, 0, y) + G(0, y)H(l, 1, 0, z) 

+ | (G(l - z, 1, 0, 0, y) - (1, 1,0, 1, z) - tf(l, 1, 0, z) + F(l, 0, 1, 0, z) 

- H(l, 0, 0, 1, z) - H(l, 0, z)G(l - z, 0, y) - G(0, -z, 1 - z, 0, y) 

- G(0, -z, 0, 1 - z, y) + G(0, -z, 0, y)H(l, z) + F(0, 0, 1, 0, zj) 

+ 2 (G(-z, 1 - z, 1 - z, 0, y) + G(-z, 1 - z, 0, 1 - z, y) 
+ G(-z, 0, 1 - z, 1 - z, y) + G(l - z, 1 - z, -z, 1 - z, y) - G(l - z, 1 - z, 1, 0, y) 
+ H(l, 1, 1, 0, z) + H(l, 1, 0, z)G(-z, y) + 1, z)G(-z, 0, y) 

- 0, 1, z)G(-z, y) - 0, z)G(-z, 1 - z, y) + (1, 0, z)G(l - z, 1 - z, y) 

- H(l, z)G(-z, 1 - z,0,y) + H(l, z)G(-z, -z,y) - fT(l, z)G(l - z, 1 - z, -z,y) 

- z)G(-z, 0, 1 - z, y) + G(0, -z, 1 - z, 1 - z, y) - G(0, -z, 1 - z, y)H(l, z) 
+ G(0, -z, y)H(l, 1, z) - iJ(0, 1, 1, z)G(-z, y) + H(0, 1, z)G(-z, 1 - z, y) 

- iJ(0, 1, z)G(l - z, 1 - z, y) + H(0, z)G(-z, 1 - z, 1 - y)) 

+ 3(i7(0, 1, z)G(l - z, -z, y) + (0, 0, 1, z)G(l - z, y)) 

+ 4^G(-z, -z, -z, 1 - z, y) - G(-z, -z, 1 - z, 1 - z, y) 

- G(-z, 1 - z, -z, 1 - z, y) - G(l - z, -z, -z, 1 - z, y) - G(l - z, 1 - z, 0, y) 

- G(l - z, 0, 1 - z, y) - H(l, 1, z)G(-z, -z, y) - H(l, z)G(-z, -z, -z, y) 



- 29 - 



+ H(l, z)G(-z, -z, l-z,y) + H(l, z)G(-z, 1 - z, -z, y) 

+ H(l, z)G(l - z, -z, -z, y) + H(l, z)G(l - z, 0, y) - G(0, l-z,l-z,y) 

+ G(0, 1 - z, y)H(l, z) - H(0, 0, 1, z)G(-z, y) + H(0, 1, 1, z) 

- H(0, 1, z)G(-z, -z, y) - H(0, z)G{\ - z, 1 - z, y) - G(0, y)# (1, 1, z)) 

+ 6( - G(l - z, 0, 0, y) + fT(0, 0, 1, zj) + s(g(-z, l-z,l-z,y) 
+ H(l, 1, z)G(-z, y) - H(l, z)G(-z, 1 - z, y)) 



+ 



36 



: ;() ( tmmM + G( _,, 1 - 0, y) + G(-z, 0, 1 - z, y) 



- (1, z)G(-z, 0, y)J +39(- G(-z, -z, 1 - z, y) + G(l - z, 1, 0, y) 
+ G(l - z, 0, 0, y) - H(l, 0, 0, z) + G(0, 1 - z, 0, y) - iT(0, 1, 0, z) 

- (0, 0, 1, z) + (0, 0, z)G(l -z,y)- G(0, 0, y)H(l, z)) 

+ 63 ( - H(l, 1, 0, z)G(l - z, y) + G(l, 0, 0, y) - G(0, 0, 1 - z, y)) 

+ 33(i(fl"(0, z)G(-z, 1 - z, y) - H(0, z)G(l - z, 0, y) - (1, 0, z)G(-z, y)) 

+ G(l, 1, 0, y) + 0, z)G(l - z, y) - 2fT (0, 1, z)G(0, y) 
G(0,y)F(l,0,z 



+ 



6if(0 , M) ) + 3 ( + ^W£iM^j + 2)G( _ Z , ,)) 

_ m ( _ ff (1 , , M) + ff (o, i, ,) G( i - ,)) + 

+ — ( - G(l - z, 1 - z, y) - H(l, 1, z) + z)G(l - z, y)J 

+ 102 ( - G(0, 1 - z, 0, y) + G(0, 0, y)il(l, z)) + 45(^(G(0, 1 - z, y)#(0, z) 

+ (0, 1, z)G(-z, y)) - i?(0, 1, 0, 1, z) - (0, 0, 0, 1, z)) + 

87 / \ 
+ 48G(0, 1, 0, y) + — (G(0, -z, 1 - z, y) - G(0, -z, y)tf (1, z)J 

- ^ z)G(l - z, -z, y) - G(l - z, -z, 1 - z, y) - G(0, y)) 

+ 133(G(0, 1 - z, y) + G(l - z, 0, y) - G(0, y)tf (1, z)) 
12Lff(l,0,z) 151tf(0,0,z) 177#(0, 0, z)G(l - z, y) 



+ 



4 2 
635G(l,0,y) 127G(0, y)H(0, z) 



+ 115#(0,z)G(l-z,y) + 



-80G(0,0,y) 
459F(0,z) 



313( - z)G(-z, y) + G(-z, 1 - z, y)) 

G(l - z, 0, y) - G(l, 0, y) + F(0, 1, z)) + G(l - z, y) 



+ 



3G(l-z,l,y) 21 



y ( - H(l, 1, z) - ff(0, 1, z)) + 3( - G(l, 1, y) + G(0, 1, y)) 



- 30 - 



- 6H(0, z)G(l -z,y) + 42G(0, y)H(0, z) - 18G(0, y)H{l, z) 

- y (tf (1, z)G(l -z,y) + 2G(-z, l-z,y)- 2H(1, z)G(-z, y)) - 11G ^ y) 
+ 12G(1 - z, 1 - z, y) + 66 ( - G(0, 0, y) - #(0, 0, z)) + 39 ( - (0, z)G(0, y) 



G(l-z,y) + H(l,z) + 



+ | (#(1, 0, z) - G(l, 0, y)) - | (G(0, y) + (0, z 



G(0, l-z,y)\ 135 



)-— (-tf(0,0,z)-G(0,0,y)) 



+ ± - 2) - 2G(-z, -z, 1 - z, y) + G(l - z, -z, 1 - z, y) 

- H(l, 0, z)G(-z, y) + 2H(1, z)G(-z, -z, y) - (1, z)G(l - z, -z, y) 

- G(0, -z, 1 - z, y) + G(0, -z, y)il(l, z) + G(0, 1 - z, y)iJ(0, z) 

+ fT(0, 1, z)G(-z, y) - #(0, 1, z)G(l - z, y) + fT(0, z)G(-z, 1 - z, y) 

- J ff(0,0,z)G(l-z,y) 

16G(0,y) 



+ 



8yV ; 



if (1, 0, z) - 3Jf (0, 0, z) - 2G(0, 0, y) 



+ 



+ 



10ztf(0,z) 37 

y 

55#(0, z) 



- y (H(0, z)G(1 - z, y) + tf(l, z)G(-z, y) - G(-z, 1 - z, y)) 
+ 2(G(0,0,y)+i?(0,0,z)) - 2^ 2 



5(l-y)G(0,y) 1 
8z 8y 



+ — (2z-l-2y 



G(l - z,-z,l - z,y) 



+ H(l, z)G(l - z, -z, y) + G(0, -z, 1 - z, y) - G(0, -z, y)F(l, z) 
- G(0, 1 - z, y)tf(0, z) + H(0, 1, z)G(l - z, y) + fT(0, 0, z)G(l - z, y) 
1/9,1 „ N 9 , 5 v — 



+ 



( y A n \ 9 / 



7yz 5y z 
--4. + 5y + -f--^ 



) - 24vr 2 ) x 



36V2y v z ' z 2K 2 " 2 2 

G(l - z, 0, y) - G(l, 0, y) + G(0, 1 - z, y) + H(Q, 1, z) - G(0, y)tf (1, z) 

+ #(l,z)G(-z,y) -G(-z,l-z,y) 
1 



+ 



8y 



4 

+ ft ( 5 " Wz + 5z2 ) + 4( - 4 - 6z + lOz 2 )) 
36 V y 6 y z J 



G(l, 0, y) - H(0, 1, z) - G(0, y)H(0, z) + - ( 39iT(0, z)G(l - z, y) 

)" 



237F(0, 4 + 69ff(l, z)G{-z,y) _ 69G(-z, 1 - z, y) 



H(0,z)G(l-z,y) 



, \ H ( l i z )y , G(l-z,y)y ' 
+ H(l, z)G(-z, y) 1 G(-z, 1 - z, y) 



. , - , 2 , 45 48 

+ — — — + — (1 - z 2 ) + — 



1 /9965 3 
36V 24 + 2yz' 



60 85tt 2 _ 45y\ 



- 31 - 



G(l-z,y)-H(l,z) 



2 

+ ^ (l - 2z + z 2 + 2y - 2yz) G(l - z, y) + 



19ff(0,z) 

48y 2 V~ ~" 1 " 1 ' 48(1 - z) 

1 / 260437 1 , 269 9 269z 39ztt 2 , 48 



99C 4 319C 3 



439tH 
2 



8 / 
55 / 

— ( + H(0, z)G(l -z,y) + H(0, 1, z) + 2tf (1, z)G(-z, y) 



- H(l, z)G(0, y) + G(l - z, 0, y) - 2G(-z, 1 - z, y) + G(0, 1 - z, y) - G(l, 0, y)) 
+ ^ ( - (1, z) + G(l - z, y)) + ^ (i?(0, z) + G(0, y)) + ^ - 2( 3 



16499 55(1 - z) 
864 + 



48y 



(B.13) 



7 



\ ( - G(l, 1 - z, -z, 1 - z, y) + G(l, 1 - z, -z, y)tf (1, z) + H(l, 1, 0, 1, z) 

+ 1, 0, z)G(l, y) + G(l, 0, -z, y)H(l, z) - G(l, 0, -z, 1 - z, y) 

- 0, 1, z)G(l, y) - H(l, 0, 0, 1, z) - (1, 0, z)G(l, 1 - z, y) 

- G(l, 0, y)H(l, 0, z) + G(0, -z, 1 - z, 0, y) - G(0, -z, 1 - z, y)tf (0, z) 
+ G(0, -z, 0, 1 - z, y) - G(0, -z, 0, y)tf (1, z) + G(0, -z, y)il(l, 0, z) 

+ G(0, 1 - z, -z, 1 - z, y) - G(0, 1 - z, -z, y)tf(l, z) - G(0, 1 - z, 1, 0, y) 

- G(0, 1 - z, 0, y)tf (0, z) - G(0, 1,1 -z, 0, y) - G(0, 1, 0, 1 - z, y) 
+ (0, 1, 0, z)G(l,y) + G(0, 1, 0, y)il(l, z) - G(0, 1, 0, y)tf (0, z) 

+ fT(0, 1, z)G(l, 1 - z, y) - H(0, 1, z)G(l, 0, y) + fT(0, 1, z)G(0, -z, y) 
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